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The CCFR collaboration reports a precise measurement of electroweak parameters derived from 
the ratio of neutral-current to charged-current cross-sections in neutrino-nucleon scattering at the 
Fermilab Tevatron. This ratio of cross-sections measures the neutral current couplings to quarks, 
which implies a determination of sin^ ^^n-sheii) _ q 2236±0.0028(expt.) ±0.0030(model) for rritop = 
175 GeV, ruHiggs = 150 GeV. This is equivalent to Mw = 80.35 db 0.21 GeV. The good agreement 
of this measurement with Standard Model expectations implies the exclusion of additional wqq 
contact interactions at 95% confidence at a mass scale of 1-8 TeV, depending on the form of the 
contact interaction. 


PACS numbers: 13.15.Jr, 12.15.Mm, 14.80.Er 
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In the early 1980’s, accurate measurements of neutrino 
neutral-current scattering cross-sections provided key in¬ 
put to the Standard Model’s predictions of the W and Z 
boson masses. Even with the production of copious on- 
shell W and Z bosons at high luminosity pp and e“'"e“ 
colliders, neutrino-nucleon (lyN) scattering still provides 
a measurement of electroweak parameters, in particular 
Mw/Mz, with comparable precision. More importantly, 
the high precision comparison among these distinct elec¬ 
troweak processes differing in by more than two orders 
of magnitude provides a critical test of the theory and 
the possibility to search for non-Standard Model contri¬ 
butions with very high mass scales or low probabilities 
[0,||. The measurement presented here represents the 
most precise determination of sin^ 9]^ from i/N scattering 
to date and supersedes the previous result from CCFR 
due to increased statistics and improved evaluation of 
systematic errors. 

The neutral-current (NC) and charged-current (CC) 
I'M deep inelastic scattering differential cross-sections on 
an isoscalar target of light quarks are related by 



where ul,r and dL,R are the left(L) and right(R)-handed 
couplings of the to up and down quarks, respectively. 
Small corrections to Eqn. arise from massive quark pro¬ 
duction suppression, CKM matrix effects, higher-twist 
processes, electromagnetic and electroweak radiative cor¬ 
rections, and from any isovector component in the tar¬ 
get, including heavy quark seas. Within the Standard 
Model, these left and right-handed couplings are given 
t’y -^Weak “ Qem sin^ 9^ and -Qem sin^ 0^, respectively, 
allowing a measurement of sin^ 9w from ratios of NC to 
CC, and v to V CC cross-sections. Furthermore, if the 
expression for in Eqn. H is used to extract sin^ 9w, 
it is almost equal to sin^ 9w in the “on-shell” renormal¬ 
ization scheme (sin^ ^(Mi-sheii) = to all or¬ 

ders), independent of Wtop and mniggs Uo- Therefore, 
the measurement of sin^ 9\y from vN scattering, com¬ 
bined with the precise measurements of Mz from LEP 
[^, implies a measurement of Mw- In addition, direct 
extraction of ul,r and dL,R also allows a search for a 
variety of non-Standard Model processes through com¬ 
parison of the measurements from lyN scattering with 
Standard Model expectations. 

The CCFR detector consists of an 18 m long, 690 ton 
neutrino target calorimeter with a mean density of 
4.2 g/cm^, followed by an iron toroid spectrometer. 
The target calorimeter consists of 168 iron plates, 
3m X 3m x 5.1cm apiece. The active elements are liquid 
scintillation counters spaced every two plates (10.2 cm 
of steel) and drift chambers spaced every four plates 



FIG. 1. Event length in data and Monte Carlo. The predic¬ 
tion for NC, CC i/^(77^), and CC v'ei'i'e) interactions is shown. 
Inset are comparisons in the region of the length cut of data 
and Monte Carlo with systematic errors shown. 

(20.4 cm of steel). There are a total of 84 scintillation 
counters and 42 drift chambers in the target. The neu¬ 
trino target is approximately isoscalar, with a 5.67% neu¬ 
tron excess. The toroid spectrometer is not directly used 
in this analysis. 

The Fermilab Tevatron Quadrupole Triplet neutrino 
beam is created by decays of pions and kaons produced 
when 800 GeV protons hit a BeO target. A wide band 
of secondary energies is accepted by downstream focusing 
magnets. The target is located about 1.4 km upstream of 
the neutrino detector. The target and focusing train are 
followed by a 0.5 km decay region. The interactions of the 
beam are predominantly from muon neutrinos (86%) and 
anti-neutrinos (12%), but also include a small fraction 
of electron neutrinos (2.3%). The mean energies of the 
and r'e(Pe) events are 165, 135 and 160 GeV, 
respectively. The mean q'^ exchanged in the neutrino 
interactions used in this analysis was —35 GeV^. 

Neutrinos are observed via their NC and CC interac¬ 
tions, both of which are selected in this analysis from the 
energy transferred to the struck quark which appears as a 
hadronic shower in the target calorimeter. The hadronic 
energy is measured by the variable Acai, which is the sum 
of energies observed in the first 20 scintillation counters 
(2.1m equivalent of steel) in the target downstream of the 
interaction vertex. NC events usually have no final-state 
muon and deposit energy over a range of counters typical 
of a hadronic shower (5 to 20 counters), CC events 
are distinguished by the presence of a muon in the final 
state which deposits energy typical of a minimum ioniz¬ 
ing particle in a large number of consecutive scintillation 
counters downstream of the hadronic shower. 

A “length” is defined for each event as the number 
of counters between the interaction vertex and the last 
counter consistent with the energy deposition expected 
from a single muon passing through the calorimeter. 
Events with a “short” length are identified as NC candi¬ 
dates. The separation between short and “long” events 
is made at 20 counters (2.1m of steel) for events with 
Acai < 45 GeV, 25 counters (2.6m of steel) for events 
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with 45 < i?cai < 100 GeV, and 30 counters (3.1m of 
steel) for events with EcaX > 100 GeV. As is shown in the 
length distributions in Figure |^, NC interactions lie in a 
clear peak, well-below the length cut, with a continuous 
band of GG interactions under the peak. The CC 
background is calculated to be 10.5%, 21.3%, and 21.2% 
of the “short” NC candidates in the three Eca\ regions, 
respectively. 

The data used in this analysis were taken between 1984 
and 1988 in FNAL experiments E-744 and E-770. Events 
were required to have EcaX > 20 GeV to ensure full effi¬ 
ciency of the trigger |^. Fiducial cuts were made on the 
location of the neutrino interaction in the calorimeter to 
ensure that events were neutrino-induced, that a separa¬ 
tion between long and short events could be made, and 
that events originated in the central part of the calorime¬ 
ter to maximize containment of wide-angle muons and to 
minimize the ratio of electron to muon neutrinos. The 
resulting data sample consisted of 8.1 x 10® events, and 
from these events the ratio 


R 


meas 


# of short events 

# of long events 


0.4151 ±0.0010 


was calculated. 

A detailed Monte Carlo was used to determine elec- 
troweak parameters from i?meas- The only undetermined 
inputs to this Monte Carlo were the neutral current quark 
couplings which were then varied until the Monte Carlo 
predicted an i?meas which agreed with that observed in 
the data. For the extraction of sin^^vV) the couplings 
in the Monte Carlo were fixed to their Standard Model 
predictions as functions of sin^ 9w, which was then var¬ 
ied as the only free parameter. The Monte Carlo in¬ 
cluded detector response and beam simulations, as well 
as a detailed cross-section model which included electro¬ 
magnetic radiative corrections, isovector target correc¬ 
tions, heavy quark production and seas, the longitudinal 
cross-section, and lepton mass effects. 

There are three major uncertainties in the comparison 
of i?meas from the Monte Carlo to the data: the statistical 
error in the data, the uncertainty in the effective charm 
quark mass for charged current charm production, and 
the uncertainty in the incident flux of i^e’s on the detec¬ 
tor. 

The charm quark mass error comes from the uncer¬ 
tainty in modeling the mass threshold of the charm pro¬ 
duction cross section. The Monte Carlo uses a slow- 
rescaling model with the parameters extracted using 
events with two oppositely charged muons (e.g., vq 
fi~c, c from this experiment |Q|. The shape and 

magnitude of the strange sea were also extracted in the 
same analysis and were used in the Monte Carlo cross- 
section model. This error dominates the calculation of 
Rmeas at low E^ (and low Eced) where the threshold sup¬ 
pression is greatest. 


SOURCE OE UNCERTAINTY 

5 sin^ 6w 

data statistics 

0.0019 

Monte Carlo statistics 

0.0004 

TOTAL STATISTICS 

0.0019 

flux 

0.0015 

Transverse Vertex 

0.0004 

Energy Measurement 


Hadron Energy Scale (1%) 

0.0004 

Muon Energy Loss in Shower 

0.0003 

Muon Energy Scale (1%) 

0.0002 

Event Length 


Hadron Shower Length 

0.0007 

Counter Efficiency and Noise 

0.0006 

Vertex Determination 

0.0003 

TOTAL EXP. SYST. 

0.0019 

Charm Production, s 


(me = 1.31 ±0.24 GeV) 

0.0027 

Higher Twist 

0.0010 

Longitudinal Cross-Section 

0.0008 

Charm Sea, (±100%) 

0.0006 

Non-lsoscalar Target 

0.0004 

Structure Functions 

0.0002 

Rad. Corrections 

0.0001 

TOTAL PHYSICS MODEL 

0.0030 

TOTAL UNCERTAINTY 

0.0041 


TABLE I. Uncertainties in the extraction of sin^ 9^ shell) 
from the CCER data 


The h'e(7'e) flux uncertainty has an important effect on 
.Rmeas bccause almost all charged current zze(Ve) events 
are short events. Therefore, the relatively small frac¬ 
tional uncertainty in the Ueii/e) flux is a significant effect, 
particularly at high Eca\ since most Ve{pe) charged cur¬ 
rent interactions deposit the full incident neutrino energy 
into the calorimeter. Two techniques were used to deter¬ 
mine the r'e(Ue) flux. In both E744 and E770, a detailed 
beam Monte Carlo was used to predict the flux, up to a 
4.1% uncertainty in each experiment. This 4.1% is domi¬ 
nated by a 20% production uncertainty in the K^ content 
of the secondary beam which produces 16% of the t'e flux. 
The bulk of the flux comes from decays, which are 
well-constrained by the observed spectrum from K ^2 
decays [|j. In E770, the Vf.{ve) flux was measured di¬ 
rectly using the fact that CC Ve{ve) interactions will have 
a high fraction of their energy deposited in the first three 
counters downstream of the event vertex. This gave an 
independent measurement of the Ueive) flux with a un¬ 
certainty of 3.5% which was in good agreement with the 
Monte Carlo method 1^. Combining these techniques, a 
measurement of the Ve{pe) flux sum in E744 and E770 is 
obtained with a 2.9% error. 

Other sources of experimental uncertainties were kept 
small through extensive modeling based on neutrino and 
testbeam data ^-|l^. The cross-section model used 
a modified Buras-Gaemers parameterization II of the 
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cross-hatched band. CDF measurements provide the hollow 
ellipse, and the DO measurements are shown in the striped 
ellipse. 


CCFR data for input parton distributions. This resulted 
in partial cancellations of certain systematic effects, such 
as errors in energy calibration. Systematic uncertainties 
associated with the measurement of i?cai include pos¬ 
sible small NC/CC shower differences (constrained by 
a LEPTO Monte Carlo study 0), uncertainties in the 
muon energy deposit within the hadron shower, uncer¬ 
tainties in the resolution function, c/tt response, and ab¬ 
solute energy scales obtained from hadron and electron 
test beam measurements |^,|^ . The length uncertainties 
include those associated with the shower length parame- 
terizations of test beam measurements lll[| , the calorime¬ 
ter longitudinal vertex determination (studied using the 
vertex from events with two muon tracks), counter ineffi¬ 
ciencies and noise, and counter spatial dimensions. In the 
cross-section model, the level of the charm sea was taken 
from the CTEQ4L parton distribution functions and was 
assigned a 100% uncertainty. Our parameterization of 
^long = O'/,/or is based on QCD predictions and data 
from charged lepton scattering experiments |p^ and is 
varied by 15% of itself to estimate uncertainties. A cor¬ 
rection for the difference between u and d valence quark 
distributions in nucleons, obtained from muon scatter¬ 
ing data Q, was applied to account for the 5.67% ex¬ 
cess of neutrons over protons in the target. A correc¬ 
tion was also applied for the asymmetry in the u and d 
sea distributions suggested by the NA51 Drell-Yan data 
and the Gottfried Sum Rule as measured in muon 
scattering |j^. Radiative corrections to the scattering 
cross-sections were applied using computer code supplied 
by Bardin 0 and 1-loop electroweak radiative correc¬ 
tions as calculated by Marciano and Sirlin Q. Possible 
higher-twist corrections were parameterized in a VDM- 
based model of Pumplin which was constrained by 
lepto-production data. Table || shows the uncertainties 
in the determination of sin^ 6\y. 

The extraction of sin^ 9^^ by comparing i?meas 

in the data to the Monte Carlo with sin^ 0^ as the 
single free parameter yields, 

sin^ ^ Q 2236 ± 0.0028((expt.) ± 0.0030(model) 



FIG. 3. One-sigma constraints on the isoscalar neutral cur¬ 
rent quark couplings, and gji, from this result (hatched) 
and from other neutrino data (solid). 


Interaction 

A+ 

A" 

LL 

4.7 TeV 

5.1 TeV 

LR 

4.2 TeV 

4.4 TeV 

RL 

1.3 TeV 

1.8 TeV 

RR 

3.9 TeV 

5.2 TeV 

VV 

8.0 TeV 

8.3 TeV 

AA 

3.7 TeV 

5.9 TeV 


TABLE II. 95% Confidence Lower Limits on mass scales 
of new lyi/qq contact terms from CCFR. 


-hO.0006 X 


/(m?„p-(175 GeV)2)\ 
(100 GeV)2 j 


-0.0002 X 



^Higgs \ 

150 GeV/ ■ 


( 2 ) 


The explicit dependence of the central value of this re¬ 
sult on rric is 0.2236-1- 0.0111 x (toc — 1.31 GeV). Only 
data with Ecai > 30 GeV was used in this result to 
reduce the effect of the charm-production and higher- 
twist systematics which are largest at low Ecax- Since 
sin^ = 1 — this result implies Mw = 

80.35 ± 0.21 GeV. This value of Myy agrees with direct 
mass measurements as shown in Figure y, and this result 
is also in good agreement both with previous vN mea¬ 
surements and with Standard Model expectations |^. 

Rmeas is also uscd to extract a constraint on the cou¬ 
plings of quarks to the Z'^: 


K = 0.5820 ± 0.0041 = 1.789751 + l-14795|j 

-0.0916di - 0.0782(5|j, (3) 

where 52 ^^ = + dl and Slj^ = The 

explicit dependence of the central value of k on rric is 
K = 0.5820 — (rric — 1.31) * 0.0111 The Standard Model 
prediction is k = 0.5817 ± 0.0013 for the measured values 
of mz, mtop, rnw- Figure 3 compares this result to other 
neutrino data m and the Standard Model prediction [|. 

Because the GGFR result for neutral-current quark 
couplings are in good agreement with Standard Model 
expectations, this result disfavors the introduction of ad¬ 
ditional processes contributing to the same final state. To 
quantify this, we use the GGFR data to constrain models 
of new viyqq contact interactions. Such models are often 
used to parameterize searches for fermion compositeness 
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or non-zero size ]^ . We assume a generation-universal 
interaction of the form 

dz47r - u 7 — 

Such interactions would shift one or more of the quark 
couplings in 1 / or 17 induced interactions from its Standard 
Model value Q. One-sided 95% confidence level lower 
limits for each A are set by finding the points in the space 
of measured couplings at which the is 1-64 units above 
the minimum, and then determining the A to which these 
points correspond in each model. Limits for LL, LR, RL, 
RR, VV (vector, All = ^lr = ^rl = ^rr) and AA 
(axial-vector. All = -Alr = -A^l = Arh) are shown 
in Table These limits are roughly comparable to the 
limits for charged lepton-quark interactions in pp collider 
data from the Fermilab Tevatron . 

In summary, CCFR has produced the most precise 
measurements of neutral-current neutrino-nucleon inter¬ 
actions to date, and has used these measurements to 
constrain neutral-current coupling to quarks. Within 
the Standard Model, this leads to a measurement of 
sin^ ^ ^ ^236 ± 0.0041 (mtop = 175 GeV, 

WHiggs = 150 GeV) which corresponds to Mw = 80.35 ± 
0.21 GeV. This result is also used to limit possible TeV- 
scale contact interactions of neutrinos and quarks outside 
the Standard Model. 
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